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3. (Amended) Method for the simplification of a source mesh according to claim 
1 , characterized in that it comprises a step for the selection of an edge merger to be made 
among all the edge mergers possible, taking account of: 

- at least one piece of information representing the curvature defined locally 
around the edge considered ; 

- at least one piece of information representing the geometrical dynamics defined 
locally. 

6. (Amended) Method for the simplification of a source mesh according to claim 
1, characterized in that said information representing the geometrical dynamics belongs to 
the group comprising : 

. - the length of the edge considered ; 

- a mean of the surfaces of the faces neighboring said edge considered ; 

- a mean of the lengths of the edges adjacent to the vertices forming said edge 
considered ; 

- a combination of the lengths of edges and/or surfaces of faces ; 

7. (Amended) Method for the simplification of a source mesh according to claim 
1, characterized in that the decimation is interrupted as a function of one of the criteria 
belonging to the group comprising : 

- a compression rate achieved ; 

- a geometrical complexity achieved, expressed by a number of vertices or faces ; 

- a threshold of curvature achieved. 
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8. (Amended) Method for the simplification of a source mesh according to claim 
1, characterized in that it constitutes a step of initialization of a method of geometrical 
optimization of a mesh. 

9. (Amended) Method for the geometrical optimization of a source mesh 
comprising a step of initialization implementing the method of simplification according 
to claim 1 . 

10. (Amended) Method for the encoding of a source mesh (M) according to claim 
9, representing a 3D object, delivering a simplified mesh (M') corresponding to said 
source mesh (M), said meshes being defined by a set of vertices, edges and/or faces, 
characterized in that it implements the method for the simplification of a source mesh (M) 
representing a plurality of surfaces defined by vertices, faces and orientations of these 
faces, said method implementing a step of decimation by edge merger, consisting of the 
association of an edge to be decimated, defined by two vertices, with a single vertex so as 
to obtain a simplified mesh M', 

characterized in that the method comprises a pseudo-optimizing step after said 
step of decimation by merger of a edge, positioning the vertex resulting from said merger 
so as to reduce the geometrical deviation between said source mesh M and said simplified 
mesh M', and then a step of minimization of a volume contained between said source 
mesh (M) and said simplified mesh (M'). 

12. (Amended) Method for the encoding of a source mesh according to claim 10, 
characterized in that said minimizing step implements an iterative process progressively 
optimizing the positions of the vertices of said simplified mesh (M'). 
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14. (Amended) Method for the encoding of a source mesh according to claim 1 1, 
characterized in that said step of minimization implements an adaptive gradient method. 

16. (Amended) Method for the encoding of a source mesh according to claim 10 
characterized in that, at each iteration, an elementary variation of said volume 
corresponding to a vector field SM is determined and in that, since the surface is 
parametrized by u and v so that a vector SM is expressed in the form SM (u, v), said 
elementary variation is likened to the parallelepiped generated by the evolution of the 
surface area element dudv in the direction SM (u, v). 

17. (Amended) Method for the encoding of a source mesh according to claim 10, 
characterized in that said simplified mesh is parametrized by means of a model of finite 
elements. 

19. (Amended) Method for the encoding of a source mesh according to claim 10, 
characterized in that said method implements a progressive encoding of said simplified 
mesh by decimation and local optimization. 

20. (Amended) Method for the encoding of a source mesh according to claim 10, 
characterized in that it comprises a step of limitation of the deterioration due to an 
elementary conversion implementing a priority queue on the elementary conversions. 

23. (Amended) Application of the method for the encoding of a source mesh 
according to claim 1 to at least to the following fields : 

- virtual reality ; 

- scientific simulation ; 

- modelling. 
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REMARKS 



The above preliminary amendment is made to remove multiple dependencies from 
claims 3, 6, 7, 8, 9, 10, 12, 14, 16, 17, 19, 20 and 23. 

Applicants respectfully request that the preliminary amendment described herein 
be entered into the record prior to calculation of the filing fee and prior to examination and 
consideration of the above-identified application. 



communication, please contact Applicants' primary attorney-of record, John J. Gresens (Reg. No. 
33,112), at (612) 371.5265. 



If a telephone conference would be helpful in resolving any issues concerning this 



Respectfully submitted, 



MERCHANT & GOULD P.C. 
P.O. Box 2903 

Minneapolis, Minnesota 55402-0903 
(612) 332-5300 



Dated: March 28, 2001 
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Reg. No. 33,112 
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3. (Amended) Method for the simplification of a source mesh according to [any of 
the claims 1 and 2] claim 1 , characterized in that it comprises a step for the selection of 
an edge merger to be made among all the edge mergers possible, taking account of : 

- at least one piece of information representing the curvature defined locally 
around the edge considered ; 

- at least one piece of information representing the geometrical dynamics defined 
locally. 

6. (Amended) Method for the simplification of a source mesh according to [any of 
the claims 1 to 5] claim 1 , characterized in that said information representing the 
geometrical dynamics belongs to the group comprising : 

- the length of the edge considered ; 

- a mean of the surfaces of the faces neighboring said edge considered ; 

- a mean of the lengths of the edges adjacent to the vertices forming said edge 
considered ; 

- a combination of the lengths of edges and/or surfaces of faces. 

7. (Amended) Method for the simplification of a source mesh according to [any of 
the claims 1 to 6] claim 1 , characterized in that the decimation is interrupted as a function 
of one of the criteria belonging to the group comprising : 

- a compression rate achieved ; 

- a geometrical complexity achieved, expressed by a number of vertices or faces ; 



- a threshold of curvature achieved. 

8. (Amended) Method for the simplification of a source mesh according to [any of 
the claims 1 to 7] claim 1, characterized in that it constitutes a step of initialization of a 
method of geometrical optimization of a mesh. 

9. (Amended) Method for the geometrical optimization of a source mesh 
comprising a step of initialization implementing the method of simplification according 
to [any of the claims 1 to 7] claim 1 . 

10. (Amended) Method for the encoding of a source mesh (M) according to claim 
9, representing a 3D object, delivering a simplified mesh (M') corresponding to said 
source mesh (M), said meshes being defined by a set of vertices, edges and/or faces, 
characterized in that it implements [a step of simplification according to any of the claims 
1 to 7] the method for the simplification of a source mesh (M) representing a plurality of 
surfaces defined by vertices, faces and orientations of these faces, said method 
implementing a step of decimation by edge merger, consisting of the association of an 
edge to be decimated, defined by two vertices, with a single vertex so as to obtain a 
simplified mesh M', 

characterized in that the method comprises a pseudo-optimizing step after said 
step of decimation by merger of a edge, positioning the vertex resulting from said merger 
so as to reduce the geometrical deviation between said source mesh M and said simplified 



mesh M', and then a step of minimization of a volume contained between said source 
mesh (M) and said simplified mesh (M'). 

12. (Amended) Method for the encoding of a source mesh according to [any of 
the claims 10 and 11] claim 10 , characterized in that said minimizing step implements an 
iterative process progressively optimizing the positions of the vertices of said simplified 
mesh (M'). 

14. (Amended) Method for the encoding of a source mesh according to [any of 
the claims 1 1 to 13] claim 11 , characterized in that said step of minimization implements 
an adaptive gradient method. 

16. (Amended) Method for the encoding of a source mesh according to [any of 
the claims 10 to 15] claim 10 characterized in that, at each iteration, an elementary 
variation of said volume corresponding to a vector field SM is determined and in that, 
since the surface is parametrized by u and v so that a vector SM is expressed in the form 
SM (u, v), said elementary variation is likened to the parallelepiped generated by the 
evolution of the surface area element dudv in the direction SM (u, v). 

17. (Amended) Method for the encoding of a source mesh according to [any of 
the claims 10 to 16] claim 10 , characterized in that said simplified mesh is parametrized 
by means of a model of finite elements. 



19. (Amended) Method for the encoding of a source mesh according to [any of 
the claims 10 to 18] claim 10 , characterized in that said method implements a progressive 
encoding of said simplified mesh by decimation and local optimization. 

20. (Amended) Method for the encoding of a source mesh according to [any of 
the claims 10 to 19] claim 10 , characterized in that it comprises a step of limitation of the 
deterioration due to an elementary conversion implementing a priority queue on the 
elementary conversions. 

23. (Amended) Application of the method for the encoding of a source mesh 
according to [any of the claims 1 to 22] claim 1 to at least to the following fields : 

- virtual reality ; 

- scientific simulation ; 

- modelling. 
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METHOD FOR THE ENCODING OF A SOURCE MESH WITH 



OPTIMIZATION OF THE POSITION OF A VERTEX RESULTING FROM AN 



DUElWERG^RTAT«TX^5RRESPONDING APPLICATIONS 



1. FIELD OF THE INVENTION AND^PTO'C ATlON S 

1.1 Field of the invention 

The field of the invention is that of the encoding of geometrical data 
structures or meshes, especially large-sized meshes. More specifically, the 
invention relates to the representation and encoding of objects or scenes in 
three dimensions. More specifically again, the invention relates to a 
technique of approximation of a 3D source mesh that can be used alone or in 
combination with other known techniques. In the latter case, the method of 
the invention may be an advantageous step of initialization. 

A mesh is conventionally defined by a set of vertices and oriented 
faces defining a topology. Meshes of this kind are for example used in 
computer graphics to model 3D objects with a limited geometrical complexity. 

The approximation of a mesh M is done by finding a mesh M' that has 
a geometrical complexity that is lower than the geometrical complexity of the 
mesh M and approaches the geometry of M to the utmost possible extent. 

1.2 Exemplary applications 

The invention can be applied in all fields where it is desirable to 
reduce the number of information elements needed to represent and/or 
efficiently manipulate a 3D object or a set of objects, for example to analyze 
it, store and/or transmit it and/or provide for its rendition. 

By way of an indication, the invention can be applied especially to the 
fields of : 

virtual reality (visits or virtual shops, leisure activities, remote 
handling, etc.). In this type of application, the approximation of the 
meshes reduces the cost of rendition of complex scenes, 
especially by defining the notion of scaleability on meshes (in 
terms of the function of the viewpoint, graphic capacities, desired 
refresh rate, etc.). In the case of distributed or shared virtual 
reality, this can also be used to adapt the complexity of a scene to 
the rendition and storage capacities of the different terminals, as 
well as to the bit rates of the networks ; 
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scientific simulation (finite elements, CAO, etc.). The reduction of 
the geometrical complexity of the models brings accelerated 
computation time and faster decision-making, especially in CAO 
designing, and eliminates redundant information in a 3D database; 
5 - modelling (3D scanner with reconstruction of surfaces from non- 

organized points, volume scanners, reconstruction of surfaces 
from stereoscopic photos or video sequences, digital models of 
terrains with satellite imaging or radar, etc.). A digital terrain 
model can thus be used to obtain a mesh representing the 
10 topology of a region. A mesh of this kind is obtained by the 

regular sampling of an image storing the altitude information at 
each point. The result thereof is a large quantity of information 
comprising information that is useless for scientific simulation or 
too costly for rendition (in the case of simulators). The 



15 approximation of meshes reduces the quantity of data while 

ensuring high geometrical fidelity to the initial data and the 
preservation of the topology. 
2. PRIOR ART 
2.1 The families of algorithms 
20 Several techniques of mesh approximation are already known. The 



most widespread meshes may be classified under three great families of 
algorithms depending on whether they work by: 

decimation; 

sub-critical resampling; 
25 - adaptive subdivision. 

2.1.1 Decimation 

Decimation consists in Tteratively withdrawing vertices and/or faces 
from a mesh. This operation is called an elementary simplification operation. 
The methods that implement this principle of decimation may also optimize 
30 the positions of the vertices after or during the simplification, the latter being 
chosen so as to preserve the topology of the mesh to the utmost extent. 

2.1.2 Sub-critical re-sampling 

The re-sampling consists in sampling an original model, either by 
taking points randomly on its surface and then retriangulating or by defining a 
35 3D grid and agglomerating the vertices in each elementary box of the grid. 



The model thus generated is simplified and must approximate the initial data 
as closely as possible. This is a fast technique but preserves neither the 
topology nor the visually important characteristics of the meshes. 

2.1.3 Adaptive subdivision 

Adaptive subdivision starts with a model comprising a very simple 
geometry that is then recursively subdivided by adding a detail at each 
iteration in the regions in which the approximation error is the maximum. 

2.1 .4 Combinations of algorithms 

In order to enable an approximation of a mesh with a satisfactory 
quality of reconstruction, it is necessary to combine a decimation and an 
optimization of the positions of the vertices preserved. In other words, since 
the basic goal of a source mesh encoding method is to maximize the quality 
of the approximation for a given geometrical complexity, this goal must have 
the following properties: 
decimation; 

preservation of the topology; 
- optimization of the positions according to a predefined error 
criteria. 

There is also a first known method called "remeshing" that meets 
these criteria. It has been presented especially in the document by 
Greg Turk "Retiling Polygonal Surfaces", SIGGRAPH 92 Conference 
proceedings, pp. 55-64, 1992). It works by sampling, decimation and 
optimization of positions. It can be used to parametrize the number of 
desired levels of resolution and the final number of vertices for each of them. 

The algorithm starts with a random distribution of the' number of 
vertices parametrized on the surface of the model and then optimizes their 
position by binding them by forces of repulsion (functions of local curvature). 
It then triangulates the polygons thus formed by including the vertices of the 
original mesh. Finally, the model is decimated one vertex after the other (if 
this operation keeps the local topology) and then the resulting aperture is 
retriangulated. 

Thus, a refined model is obtained in the regions of high curvature and 
a decimated model is obtained in the plane regions. 



One drawback of this method is that it has a random part. 
Consequently, two successive decimations do not give the same result. 
Furthermore, this technique tends to smoothen the discontinuities in most 
situations. 

Furthermore, fidelity to initial data is not expressed directly during the 
optimization of the positions. 

Another technique known as "progressive mesh encoding" has been 
developed by Hugues Hoppes in "Progressive Meshes", SIGGRAPH 96 
Conference Proceedings, pp. 99-108, 1996. It relies on the decimation and 
optimization of points. This technique is described in greater detail in 
Appendix 1 with reference to Figures 24 to 27, so as not to burden the 
present discussion. It is clear however that this Appendix 1 as well as 
Appendix 2 form an integral part of the description. 

This technique of progressive mesh encoding has various drawbacks. 
In particular, it does not naturally preserve the geometrical discontinuities or 
singularities; This means that there has to be a management of special 
cases such as corners, finishing sharp edges and regular sharp edges to 
prevent any breaks in topology. This leads to complicated computations that 
are difficult to implement and assume difficult parametrizing operations (what 
is a sharp edge? how to weight the bounce term to prevent a smoothening of 
the surfaces? etc.). 

Furthermore, it is not efficient for characterizing the perceptual 
differences between two meshes. 
3. GOALS OF THE INVENTION 

The invention especially relates to the technique of decimation 
implemented for example by these different techniques of mesh encoding or 
by others. ~ - - 

The goal of the invention especially is to overcome the different 
drawbacks of the known techniques; 

Thus, a goal of the invention is to provide a method of mesh 
simplification by decimation (edge merger) that is more effective in terms of 
perceptual quality than known techniques. 

In other words, a goal of the invention is to give a method of this kind 
which, up to a high level of decimation, maintains the singularity of the 
meshes and preserves the topology. 



A goal of the invention is also to give a method of this kind that is 
simple to implement in terms of computations to be made and has high 
execution speed. 

According to a first aspect of the invention, another goal is to provide a 
method of this kind that can be used alone in order to give a fast method of 
mesh simplification. 

According to a second aspect of the invention, one goal is to give a 
method of decimation of this kind that can be used, so that is improved, in a 
method of geometrical mesh optimization. 

According to yet another aspect of the invention, a goal of the 
invention is to provide a method for the encoding of a source mesh in three 
dimensions with a better ratio between quality of approximation and 
geometrical complexity than known techniques; 

In particular, the invention is aimed at giving a method of this kind that 
more faithfully keeps the major characteristics (singularities or discontinuities) 
that have to be preserved. 

Another goal of the invention is to give a method of this kind that does 
not require searching for and managing special cases (such as edges and 
corners) or implementing specific parameters (for example criteria of 
recognition of a edge, thresholds, etc.) that have to be defined for each 
mesh. 

Yet another goal of the invention is to provide a method of this kind 
enabling a successive reconstruction of the mesh, this mesh being very 
quickly recognizable in a rough representation. 

It is also a goal of the invention to give a method by which it is 
possible to fulfil the functions laid down in the ISO-MPEG4 standardization 
project, namely: - - _ - - 

- the compression of the meshes which must greatly reduce the 
data describing the geometry, with minimum visual loses; 

- the progressive representation of a mesh, enabling the defining of 
a mesh in the form of a low-resolution basic mesh and a sequence 
of refinements; 

- the matching of a wide range of bit rates, the quality being 
acceptable, but for very low bit rates. 



4. MAIN CHARACTERISTICS OF THE INVENTION 

4.1 Study of optimization of the position of the vertices 

These goals as well as others that shall appear more clearly 
hereinafter are achieved according to the invention by means of a method for 
the simplification of a source mesh (M) representing a 3D object, delivering a 
simplified mesh (M') corresponding to said source mesh (M), said meshes 
being defined by a set of vertices, edges and/or faces, said method 
implementing a step of decimation by edge merger, wherein with a edge to 
be decimated, defined by two vertices, there is associated a single vertex so 
as to obtain a simplified mesh M'. 

This method comprises an optimizing or pseudo-optimizing step after 
said step of decimation by merger of a edge, providing for the positioning of 
the vertex resulting from said merger so as to reduce the geometrical 
deviation between said source mesh M and said simplified mesh M'. 

Thus, the particular aspects, especially the sharp edges of the mesh 
are complied with. 

Said step of pseudo-optimization may advantageously consist in 
enumerating the sharp edges around two vertices forming the edge to be 
merged and distinguishing the following two cases: 

- if the numbers of sharp edges are the same around two vertices, 
the vertex resulting from the merger is placed in the middle of the 
segment linking said vertices; 

- if the numbers of sharp edges are different, the vertex resulting 
from the merger is placed on the vertex with the greatest number 
of sharp edges. 

4.2 Advantageous characteristics of the invention 

The invention can- be "implemented especially by means of an 
algorithm organized in two parts: the optimizing of positions and decimation. 

Decimation is used to obtain a simplified geometrical mesh, from the 
original mesh, while at the same keeping the topology and a high degree of 
resemblance with this original mesh. 

The optimizing of the positions of the vertices gives the best 
approximation of the volume-based metrics of the invention. 

By combining decimation and optimization, it is possible to generate 
either a mesh in progressive form (by using alternately the decimation of a 
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vertex and a local optimization) or distinct levels of resolution (by using 
successively the decimation of a set of vertices and total optimization). 

Advantageously, the method of the invention comprises a step for the 
selection of one edge merger to be made among all the edge mergers 
possible, taking account of : 

at least one piece of information representing the curvature 

defined locally around the edge considered; 

- at least one piece of information representing the geometrical 
dynamics defined locally. 

By taking these two criteria into account it is possible, as shall be seen 
hereinafter, to optimize the choice of the mergers to be made, in eliminating 
the perceptually less significant elements as a priority. 

Preferably, said step of selection implements a queue of priorities of 
edges to be merged as a function of a priority criterion, said information 
representing the curvature and then a secondary criterion, said information 
representing the geometrical dynamics. 

This hierarchy of criteria is used to achieve high efficiency. 

Advantageously, said selection step manages a threshold of 
curvature, only the edges with a curvature smaller than said threshold being 
considered for the application of said secondary criterion, said threshold 
being increased when there is no longer any edge having a curvature below 
this threshold. 

According to other particular embodiments, said information 
representing the geometrical dynamics may belong to the group comprising: 
the length of the edge considered; 

- a mean of the surfaces of the faces neighboring said edge 
considered; ~ - - 

- a mean of the lengths of the edges adjacent to the vertices 
forming said edge considered; 

- a combination of the lengths of edges and/or surfaces of faces; 

- any other characteristic magnitude connected to the local density. 
As shall be seen hereinafter, taking into account the length of the edge 

is a simplified technique and gives very good results. 

The decimation may be interrupted in particular as a function of one of 
the criteria belonging to the group comprising : 



a compression rate achieved; 

a geometrical complexity achieved, expressed by a number of 

vertices or faces; 

a threshold of curvature achieved. 

According to a first embodiment of the invention, the method of 
simplification of a source mesh forms a step of initialization of a method of 
geometrical optimization of a mesh. 

The invention also relates to a method of geometrical optimization of 
this kind of a source mesh comprising a step of initialization implementing the 
method of simplification described here above. 

According to a second mode of implementation of the invention, the 
method of simplification of a source mesh may be used alone. 

Thus, the method of mesh encoding may implement a step of 
minimization of a volume contained between said source mesh (M) and said 
simplified mesh (M'). 

Thus, the invention relates to a technique of approximation of meshes 
using a volume-based metrics and not, as is conventionally the case, a 
metrics based on the distance between a vertex and a surface (a distance of 
this kind furthermore is not unique). 

This new approach naturally takes account of the major characteristics 
to be kept (singularities) without it being necessary to have recourse to the 
detection of special cases. No particular parametrizing is required before the 
processing of a mesh. The volume between two meshes is simply 
minimized. 

Advantageously, since each of said meshes is defined by the position 
of each of the vertices and edges that connect them, said minimizing step 
provides for the determining of the position of the vertices (X-j , X2, X n ) of 
said simplified mesh (M') minimizing the volume V(M, IvT) between said 
source mesh and said simplified mesh. 

Preferably, said minimizing step implements an iterative process 
progressively optimizing the positions of the vertices of said simplified mesh 
(M'). 

Said iterative process may especially be interrupted when at least one 
of the following stopping criteria is achieved: 

a maximum number (N) of iterations; 



a difference between two successive shift vectors of the positions 
of the vertices that is smaller than a predetermined threshold (e). 
According to an advantageous embodiment, said step of minimization 
implements an adaptive gradient method. This method, as shall be seen 
hereinafter, implements the invention with a reduced complexity of 
computation. Indeed, it is no longer necessary to carry out an explicit 
computation of volume. 

In this case, said adaptive gradient method may advantageously rely 
on the following operations: 

- the selection of a vector X p of R 3n (n > 1 ) of said simplified mesh 
and computation of the gradient VE (X p ) in X p of the function to be 
minimized E = d(M,M') : 

- the determining of the position of X p * of X p of said mesh 
according to the relationship defined in the iteration k+1 by: 

k varying from 0 to n-1 (with n < N) and Yk being the step of said 
relationship : 

N represents the number of useful iterations to attain the 
convergence (this convergence being achieved when the 
difference between two successive determinations is smaller than 
a threshold (e) or when a maximum number of iterations is 
achieved). 

It may be noted that X p belongs to R 3 if a local optimization is done or 
to R 3n if a total optimization is done. In other words, the positions of a vertex 
are optimized in the event of focal optimization and the positions of all the 
vertices are optimized simultaneously (in the form of a singe vector of R 3n ) in 
the case of a total optimization. 

Preferably, said step Yk is variable and varies as a function of the 
oscillation of two successive shifts of the vertices and/or as a function of the 
variations of the energy. It may be maintained especially between two 
terminals Ymm ar >d Ymax- 

According to an advantageous embodiment of the invention, at each 
variation, an elementary variation of said volume corresponding to a vector 
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field SM is determined. Since the surface is parametrized by u and v, so 
that a vector SM is expressed in the form 8M (u, v), said elementary 
variation is likened to the parallelepiped generated by the evolution of the 
surface area element dudv in the direction smt (u, v). 

This elementary variation of the volume between M and M' may 
advantageously be written as follows: 

dd{M,M,SIVf ) = - J J ri(u,v)n(u,v)6M(u,v)do(u,v) 

With : yj(u, v) = -n.(M, M', n ( UjV )) = 1 if the unit volume rc(u,v) is oriented 
towards the interior volume and -1 if not. 

According to a preferred embodiment of the invention, said simplified 
mesh is parametrized by means of a model of finite elements. Said finite 
elements are advantageously obtained by means of a refined interpolator. 

In this case, said simplified mesh can be written as follows: 

N 

M(u,v)=^X i ^ i (u,v) 

with X ; (u,v) being a shape function matched to a model of triangular finite 
elements; 

and Xj being the vertex of said mesh defined in R 3 , 

and a shape function is defined on the mesh by means of barycentric 
coordinates. 

By using the previous expression, we get: 

N 

od(M,M ,oX, 6X n ) = - J J rj(u,v)n(u,v) £ 6X i X 1 (u,v)da(u,v) 
= ~ L J J , n(u,v)n(u,v)6X i X. i (u,v)da(u,v) 

i = l u,v 

Advantageously, the position of a vertex Xj at the k-th iteration is 
written as follows: 

3 8d 

with the partial derivative of the distance for a vertex Xj : 
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= - J J* rj(u,v)ri(u,v)X i (u,v)do(u,v) 

dX l u ,v G Supp(X,) 

Advantageously, the method of the invention implements a 
progressive encoding of said simplified mesh by decimation and local 
5 optimization. Thus, said simplified mesh is advantageously represented by a 
basic mesh and a sequence of refinements of said basic mesh. 

Advantageously, the computation of the gradient at the k-th iteration 
comprises the following steps: 

- discretizing of the expression of the partial derivative of the 
10 distance for each vertex Xj in the form : 

^- = -ZIZ r,(i T ,j T )n(i T ,j T )A. 1 (i x ,j T )da(i T) j,) 

OJ ^i teS^ = lj t =l 

with: 

15 - S as the set of triangles neighboring the vertex Xj ; 

N x the number of points sampled in the direction of u ; 
M T the number of points sampled in the direction of v ; 
computation of the orientation of the surfaces by identification of 
the closest intersection with the source mesh M. 
20 According to a preferred aspect of the invention, the relative 

orientation of the surfaces with respect to said source mesh and said 
simplified mesh is computed according to the equation: 

ri(u,v) = -<n N , r ,n M > ,<3 M , XjX M > 

25 

with : 

Xj as the point sampled by M' ; 

Xm the point of intersection of the straight line passing through Xj and 
having the direction h with the source mesh M ; 
30 fi M the normal to the source mesh M at the point X|y] ; 

n w the normal to the source mesh M' at the point Xj ; 
^ X-X M represents a vector). 

According to another advantageous aspect of the invention, the 
method comprises a step of limitation of the deterioration due to an 
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elementary conversion implementing a priority queue on the elementary 
conversions. 

Preferably, said step of limitation of the deterioration due to an 
elementary conversion comprises the steps of: 
5 - computing a cost for each possible elementary conversion ; 

carrying out the lowest cost elementary conversion ; 
recomputing the costs of the elementary conversions modified by 
the previous elementary conversion ; 

adding the new elementary conversions created and computation 
10 of the corresponding costs. 

The cost of an elementary conversion (Tj) may be expressed by: 

C(T i (X i ,X j )) = max d,(V M ,F(X[)) 

with : 

15 Tj conversion merging two vertices Xj and Xj of said simplified 

mesh M' ; 

Xj f the vertex of said simplified mesh M' resulting from said 
conversion; 

F(Xj f ) the faces of said simplified mesh M' neighboring the vertex Xj f 
20 after said conversion; 

Vm set of the vertices of said source mesh M belonging to the faces 
having been intersected during the computation of the 
orientation of the surfaces during said minimization. 
The method of the invention can be applied to a very large number of 
25 technical fields and especially at least to the fields belonging to the group 
comprising: 

virtual reality ; ~ ~ - ' - 

scientific simulation ; 
modelling. 
30 5. LIST OF FIGURES 

Other features and advantages of the invention shall appear from the 
following description of a preferred embodiment of the invention given by way 
of simple illustrative and non-restrictive examples and the appended figures, 
of which: 

35 - Figure 1 illustrates the principle of a edge merger ; 
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- Figure 2 illustrates the principle of the priority queue combining the 
curvature and the geometrical dynamics defined locally according to the 
invention ; 

- Figure 3 illustrates the enumeration of the sharp edges around a 
5 vertex of the mesh ; 

- Figure 4 shows the principle of the pseudo-optimization between 
edge merger according to the invention ; 

- Figure 5 illustrates the principle consisting in minimizing the volume 
between two surfaces ; 

10 - Figures 6 to 11 are discussed in Appendix 2 describing the 

mathematical aspects of the invention and respectively illustrating : 

- Figure 6 : elementary variation of the mesh ; 

- Figure 7 : elementary variation of the mesh in the direction of the 

normal ; 

15 - Figure 8 : system of coordinates on a triangle ; 

- Figure 9 : shape function ; 

- Figures 10 and 11 : examples of the surface developments on 
two curvatures in two dimensions ; 

- Figure 12 illustrates the updating of the face adjacencies after 
20 merger of a edge adjacent to two faces ; 

- Figure 13 illustrates the updating of the face adjacencies of the 
merger of a edge adjacent to one face ; 

- Figure 14 illustrates the sampling on a triangle for numerical 
integration ; 

25 - Figure 15 illustrates the principle of orientation pertaining to the 

surfaces so as to prevent the formation of folds on the surface when the 
initial position is too far distant from the optimum ; 

- Figure 16 illustrates the cost of an elementary conversion ; 

- Figure 17 is a block diagram of the approximation of meshes 
30 according to the invention ; 

- Figure 18 is a block diagram of the optimization of a vertex Xj ; 

- Figures 19 to 22 illustrate the behavior of the method of the invention 
on a rectangular parallelepiped after decimation on a corner (Figure 19), a 
regular sharp edge (Figures 20 and 21) and a plane (Figure 22) ; 
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- Figure 23 illustrates the different iterations (12) performed on a 
simple example according to the invention ; 

- Figures 24 to 27 are discussed in Appendix 1 presenting the known 
technique of progressive mesh encoding, and respectively illustrate : 

5 - Figure 24 (also commented hereinafter): elementary conversion 

(edge merger) ; 

- Figure 25 : initial positions; 

- Figure 26 : updating the conversions in the vicinity affected by 
the last modification ; 

10 - Figure 27 : adding the newly created conversions. 

6. GENERAL PRINCIPLE OF THE IMPLEMENTATION OF THE 

INVENTION 
6.0 Introduction 

15 According to the main aspect of the invention, there is provided a step 

for the optimization of the position of the vertices (§ 6.2) that comes within in 
the framework of a edge merger chosen for example according to the 
technique described in §6.1. This can then be implemented within the 
framework of the methods of § 6.3 and 6.4. 

20 6.1 Geometrical simplification 

As indicated here above, the invention relates especially to a new 
technique of simplification of a 3D mesh relying on the implementation of a 
priority queue combining the local curvature and the local geometrical 
dynamics. This technique especially has the advantage of preserving the 

25 singularities on the meshes up to a high level of decimation. It furthermore 
has a valuable execution speed. 

According to one aspect of the invention, a "priority queue managing 
the edge merger topological operator is constructed. This priority queue 
combines the criteria of local curvature and of local geometrical dynamics in 

30 order make the maximum use of the degree of freedom given by the order of 
the conversions to be made on the mesh. 

The simplification of a mesh M consists of the construction of a mesh 
M' with a limited geometrical complexity that preserves a low geometrical 
deflection with M. 



15 



The algorithm of the geometrical simplification must enable the 
specifying of a geometrical resolution to the nearest vertex. For this purpose, 
an elementary topological operation of simplification is chosen with the 
following good properties : conservation of topology to a certain extent of 
5 decimation and absence of creation of holes on the surfaces and 
preservation of the orientations. 

This elementary topological operator is the edge merger as defined for 
example by Hoppe (already quoted) and illustrated in Figure 1. 

The merger of edges 10 consists in merging two adjacent vertices 1 1 

10 and 12 into one vertex 13, eliminating the two faces 14 and 15 and 
positioning the vertex 13 resulting from the merger. It will be noted that this 
conversion is reversible (with the possibility of insertion 1 6 of a vertex). 

Each elementary conversion decimates the mesh approximating M'. 
The quality of the approximation therefore deteriorates during the decimation 

15 or remains at best invariant. In order to limit the deterioration of the mesh, it 
is sought of course to carry out first of all conversions that have the least 
possible effect on the model. 

For this purpose, a priority queue is defined. This priority queue 
contains all the conversions that can be made on the mesh (namely 

20 approximately the number of edges). During the decimation, the least cost 
conversion of the priority queue is done and then eliminated from the queue. 
The cost in the neighborhood modified by the previous operation is then 
recomputed and the new potential conversions on the mesh are inserted into 
the priority queue after their cost has been computed. 

25 The geometrical singularities which are highly informative parts, must 

be preserved as long as possible during a decimation. In particular, the 
regions of high curvature appear to be highly informative. Consequently, the 
first criterion of sorting on the elementary conversions is linked to the local 
curvature around the edge to be fused. 

30 The term curvature C(Xj) around a vertex Xj is applied to the maximum 

angle between the normals to two adjacent faces around the vertex Xj. 
Then, the curvature around a edge (referenced by two vertices Xj and Xj) is 
called the mean of these criteria evaluated at each vertex. 
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A second criterion based on the length of the edge to be merged, 
which reduces the geometrical density of the mesh and gives a high 
appearance ratio on the resultant triangles, is implemented. 

It is also possible to use a habitual formula representing compactness. 
5 However, the length of the edge to be merged has the following advantages : 
giving a mesh of uniform density in regions of neighboring 
curvature ; 

preserving high compactness of the triangles since this criterion 
will tend to crate equilateral triangles ; 
10 - having a low cost of computation. 

In other words, the priority queue according to the invention is based 
on simultaneously taking the following two aspects into account : 

- a small triangle is valuable only in the highly informative region 
(namely a region of high curvature) ; 
15 - it is desirable to reduce the density of a mesh in order to reduce its 

complexity. 

The two criteria taken into account according to the invention are 
combined so as to obtain the behavior illustrated in Figure 2. This figure is a 
scale of the curvature graduated from 0 to it in radians. 
20 On the regions of low curvature 21, below the first threshold 22, the 

density is limited and the compactness obtained is reasonable since the 
edges of minimum length are merged. 

The threshold defining a low curvature is then increased (23) when 
there is no longer any conversion possible on the low curvature segment 21. 
25 Thus, on the lowest level of decimation, the curvature constraint is 

automatically relaxed in order to obtain the geometrical complexity fixed. 

The priority queue therefore has two levels of constraint organized in a 
hierarchy : the curvature thereof is the priority criterion and the density 
thereof is a secondary criterion. 

30 

6.2 Pseudo-optimization 

After an operation of edge merger, it is possible, depending on the 
essential aspect of the invention, to place the vertex resulting from the 
merger in the optimum probable position. 
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For this purpose, the notion of the sharp edge is introduced. A edge is 
sharp when the angle formed by the perpendicular to the two adjacent faces 
is greater than a fixed parametrizable threshold. The number of sharp edges 
around the vertices of the edge to be merged is then enumerated as shown 
5 in Figure 3. The vertex 31 is not associated with any sharp edge. The vertex 
32 has two sharp edges and the vertex 33 has three sharp edges. 

From this enumeration, two cases are deduced for initialization as 
shown in Figure 4 : 

if the number of sharp edges around X a and around Xb are 
10 identical, the initialization is done in the middle of the segment 

forming the edge to be merged. On the plane zones of low 
curvature, this makes it possible to preserve high compactness on 
the neighboring triangle. On a regular sharp edge (number of 
sharp edges equal to 2), this makes it possible to position the 
15 vertex close to the sharp edge which will be preserved during the 

optimization ; 

if the number of sharp edges around X a and around are 
different, the initialization is done on the vertex having the greatest 
number of sharp edges. In the most common cases, the optimum 
20 is achieved from this initial position. 

In the example of Figure 4, wherein the source mesh corresponds to a 
parallelepiped 41, it is observed that this heuristic approach places the edge : 
on the corner of the parallelepiped when the edge to be merged 
forms a corner 43 ; 

25 - on the regular sharp edge of the parallelepiped when the edge 

starts in the plane region and ends on the sharp edge 45. 
The situations where the number of sharp edges is the same around 
two vertices is illustrated at 42 and 44. 

30 6.3 Applications 

As shown here above, the technique of simplification of the invention 
can be implemented alone to offer a mesh approximation technique or as a 
step of initialization of a more complete procedure of geometrical optimization 
such as for example the one described here below. 
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6.3.1 Structure 

The invention can be implemented in a method of mesh approximation 
implementing a volume-based metric. One of the advantages of this metric is 
that it naturally takes account of the major characteristics to be preserved 
5 (the singularities) without its being necessary to have recourse to the 
detection of special cases. 

In general, the approximation of meshes generates a geometrically 
polygonal model that is simpler than the original model. The number of 
polygons necessary for the representation is thus reduced while at the same 
10 time preserving the best approximation, in the perceptual sense, of the 
original model. 

The corresponding problem may be formulated by a variational 
approach used to approximate a mesh M by a mesh M' (or &I) comprising a 
reduced number of triangles. This permits the shifting of the vertices by 
15 optimizing the position according to one or more criteria defined by means of 
an energy functional criterion to be minimized. 

The approximation of a mesh generates a model that is geometrically 
simpler than the original model. The number of polygons needed for the 
representation is thus reduced while preserving the best approximation (in 
20 the perceptual sense) of the original model. The quality of the approximation 
defines the level of perceptual resemblance for a fixed geometrical 
complexity. The invention uses a volume-based metric between the 
"approximating mesh" and the "approximated mesh". 

The principle of the algorithm is subdivided into two parts: decimation 
25 and the optimization of the positions. 

The elementary operation of decimation chosen is the edge merger, 
described here above. The optimization algorithm chosen is-an adaptive step 
gradient iterative algorithm in order to iteratively minimize the volume 
included between the model and the original mesh. 
30 By combining decimation and optimization it is possible to generate 

either a mesh progressively or distinct levels of resolution. 

6.3.2 Decimation 

In order to simplify the geometry of a mesh, it is necessary to iteratively 
decimate the mesh while at the same time keeping an appropriate topology. An 
35 elementary conversion is therefore chosen : the edge merger which merges two 
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adjacent vertices into one vertex, eliminates two faces and places the new vertex 
resulting from the merger according to the above-described principle. 

This decimation may be done alternately or successively with an optimization 
(such as the one described here above) of the positions of one or more vertices 
5 depending on the application in view. It is possible to build several levels of 
resolution of a mesh by using this elementary conversion, and the minimum grain 
sizes obtained by a decimation since the mesh is then defined to the nearest vertex. 
The method then uses an alternation of local decimation and local optimization of the 
vertex resulting from the merger. 
10 Conversely, it is possible to define levels of resolutions in a general way by 

first of all decimating a set of vertices and then comprehensively optimizing the 
positions of the vertices of the mesh. 

The successive decimation is used to construct progressive meshes. In this 
case, it is appropriate to make a careful choice of the order of performance of the 
15 conversion. This order is defined by means of a priority queue on the elementary 
conversion so as to carry out first of all the conversions that have the least possible 
effect on the original model. 

This implies defining the notion of cost for each elementary operation and the 
priority queue will be sorted out by this. Thus initially, all the possible conversions 
20 are listed, the cost of each of them is computed (without making any modifications 
on the mesh), the low-cost conversion is made and the costs of the elementary 
conversions modified by the previous conversion are recomputed. 

Naturally, techniques other than decimation may be used to obtain the 
approximate or approximated mesh, such as techniques of sub-sampling, adaptive 
25 subdivision, etc. 

6.3.3 Optimization based on volume 

The problem of optimization therefore consists in minimizing an error 
functional defined by: 

E(M, M') = d(M, M') 

30 where d(M, M') characterizes a distance between the meshes M and M' and ensures 
fidelity to the initial data. The optimization algorithm is chosen especially for its 
performance characteristics in terms of speed and/or relevance of the final solution. 

An explanation is given here below of the theoretical context of the volume- 
based metrics used to characterize the approximation error between two meshes as 
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well as the variational formulation used to optimize the quality of approximation 
according to the invention. 

6.3.3.1 Metrics 

An exemplary measurement of the error or distance between the meshes M 
5 and M' is illustrated in Figure 5. The measurement of the error is given by the 
volume V, contained between the two surfaces formed by the meshes M and M'. 
This volume is defined by the following integral of Lebesgue : 

V(M,M)= J"I v (M,M)(q)d 3 q 

s 3 

10 

with l v (q) being the function indicating the interior of the volume and q being a 
generic point of R 3 . 

6.3.3.2 Functional 

Since a mesh is partially defined by the positions of the vertices, it is possible 
15 to formulate the problem of the minimizing of the volume by means on an energy 
functional to be minimized, this functional being defined by: 

E = E error = d(M, M') = dQAM'^-XJ) 
The algorithm of optimization implemented then consists in determining the 
optimum positions of the vertices X x , X 2 , X n that minimize the distance d(M', M). 
20 6.3.3.3 Elementary variation in volume 

The elementary progress of the mesh M' is defined by the field of vectors SM . 
An expression of the elementary variation of volume is described in Appendix 2 
(§2). 

6.3.3.4 Implementation of the adaptive gradient 

25 The technique of optimization used to resolve the problem of minimization 

may advantageously rely on an iterative approach such as the' adaptive gradient. 
Appendix 2 (§ 1) gives a detailed description of the implementation of this technique 
of the adaptive gradient as well as the computation of an elementary variation of 
volume used to determine the gradient at each iteration, and the parametrizing of the 
30 mesh by a model of triangular finite elements. 

A model of triangular finite elements used to compute the gradient at each 
vertex is described in Appendix 2 (§ 3). 

The algorithm of the gradient for a vertex X ; is written in the following form : 

35 x -.xr- Yk — 
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The expression of the partial drift of the difference for a vertex Xj is: 
- J J r|(u,v)fi(u,v))v l (u,v)da(u,v) 

5 

The support of the shape function is defined by the triangles neighboring the 
vertex X;. The integral is computed by means of a discretization on the triangles 
(described hereinafter) and a heuristic method has been adopted to determine the 
relative orientation of the surfaces formed by M and M' (also described hereinafter). 
10 At each iteration of the gradient, a shift is applied to the vertex X 4 that tends 

to minimize the volume between the surfaces. The shift is the sum of the 
contributions related to the computations of orientation at several points of each 
triangle. 

The step y t is adaptive. At the outset, it is initialized at y 0 (a fraction of the 
15 length of the merged edge) and then multiplied by k (with k < 1) in the presence of 
an oscillation (namely when two vectors of successive shift have opposite directions: 

asd k+1 <0 

20 This step Vj is bounded (with a lower boundary of y min ) in order to keep the 

convergence of the algorithm. We have y £ ]8, 2&-S] ,with 0 < S -< ^ an ^ 

||HE(X)||<l/Ao (where HE represents the Hessian value of E). 

We shall consider the convergence to be achieved when the shift of the vertex 

Xj from one iteration to the other is smaller than e (e being a fraction of the local 
25 dynamics of the edges of the mesh). The number of iterations of the algorithm may 

also be limited. 

6.3.4 Description of a detailed embodiment 

Since the method of the invention can be used to generate approximations of 
meshes at different levels of resolution, a gradual encoding of meshes by decimation 
30 and local optimization may be implemented. 

To this end, an invention defines an elementary conversion of decimation 
adapted to progressive encoding, a method of local optimization and a priority queue 
on the conversions in order to maximize the ratio between quantity of information 
and visual quality at each iteration of the algorithm. Hereinafter, a description is 
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given of the implementation of this algorithm, especially the computation of the 
gradient defined on a vertex of the mesh. 
6.3.4.1 Elementary conversion 

Progressive encoding is used to define the notion of sampleability on a mesh. 
The maximum sampleability is obtained when the resolution of the mesh can be 
chosen to the nearest vertex. A mesh can thus be encoded in the form of a basic mesh 
(with a visually acceptable resolution) and a sequence of refinements. 

In order to obtain a data structure of this kind, it is necessary to interatively 
decimate the mesh, locally optimize the positions of the vertices so as to locally 
approximate the initial mesh to the best possible extent and simultaneously record 
the decimation sequence in order to encode the mesh according to the refinements. 

A reversible elementary conversion is therefore chosen, for example the 
merger of edges defined by Hoppe (document already quoted) illustrated in Figure 1. 
The edge merger 10 consists in merging the two adjacent vertices 11 and 12 into one 
vertex 13, eliminating the two faces 14 and 15 and optimizing the position of the 
vertex 13 resulting from the merger. 

This conversion is reversible (insertion 16 of a vertex). 

After conversion, the adjacencies must be updated. For this purpose, we 
consider a face/vertex type of structure and an adjacency defining all the 
relationships between the elements. The faces store the links towards their 
neighboring faces and each vertex lists the faces that contain it and the neighboring 
vertex. Figures 12 and 13 illustrate the updating of the adjacencies of the faces after 
merger of an edge adjacent to two faces (Figure 12) or one face (Figure 13). 

The decimation is thus naturally taken into account on all the zones of the 
mesh without its being necessary to detect the particular cases, such as the edges 
adjacent to a single face. 
6.3.4.2 Initialization 

The optimization acts on the position of the vertex resulting from the merger 
of the edge formed by the vertices X a and X b . In order to start the optimizing 
algorithm with a good initial condition, it is chosen to position the vertex as closely 
as possible to the probable optimum. For this purpose, the notion of a sharp edge is 
introduced. A edge is sharp when the angle formed by the normals to the two 
adjacent faces is greater than a fixed parametrizable threshold. The number of sharp 
edges around the vertices of the edge to be merged is then enumerated according to 
the approach already discussed. 
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6.3.4.3 Optimization 

It may be recalled that the expression of the functional considered: 

E = E error = d(M, M') 
This functional is minimized by optimizing the position of the vertex Xj by a 
gradient algorithm: ^ 

The determining of the step Yi and the stopping criteria are discussed in 
§ 6.3.3.4. 

6.3.4.4 Discretization 

The computation of the gradient therefore requires a discretization of the 
integral: 

= -J J ri(u,v)fi(u,v)X 1 (u,v)da(u,v) 

5X £ u,v e SuppO.) 

On the support of the shape function defined by X { . 

We are thus led to sample the surface of the triangles before assessing the 
value of the shape function Yj(x, y , z) and the orientation of the surfaces T| (u,v) at 
each point X(,x,y,z). 

The gradient is therefore expressed in the following discrete form: 

N - M 

ffr— EI L^h)mj x )h(K^o(i x ,j T ) 

with : 

S as the set of neighboring triangles of the vertex X ; ; ~ 

N T the number of points sampled in the direction u ; 

M T the number of points sampled in the direction v. 
For a given triangle, the sampling is done in the plane passing by its edges, 
the number of points being proportional to the surface of the triangle, as shown in 
Figure 14. 

The sampling reference system is formed by the edge 161 with the greatest 
length and the height 162 of the current triangle. Let n be the minimum number of 
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points per triangle fixed in advance, S min the minimum area of the triangles 
neighboring X ; and h the sampling step, h is deduced by the following formula : 

The shape function is assessed by a ratio of surfaces as described in 
5 Appendix 2. 

6.3.4.5 Orientation of the surfaces 

The relative orientation of two surfaces defined in three dimensions is 
difficult to define with exactness. Advantageously, a heuristic system is used to 
10 compute the orientation of the two surfaces M and M' (namely the term ri(u,v) of the 
equation of the gradient presented here above : 

T } (u,v) = -<n M ,n M >.<n fvf , X j X M > 

15 with: 

Xj being the point sampled on M' ; 

X M the point of intersection of the straight line passing through Xj and having 
a direction M - with the source mesh M ; 
w M being the normal of the source mesh M at the point X M ; 
20 n M ' being the normal of the source mesh M' at the point Xj ; 

where: 

<,> represent the scalar product. 

Figure 15 illustrates the computation of the orientation in two dimensions. 
In three dimensions, it is enough to replace the segments by triangles. This heuristic 
25 method greatly limits the formation of folds on the surface, since the term 

< nM,n M > expresses the opposition of the normal when they are in opposite 
directions. ~- ' 

6.3.4.6 Priorities on the conversions 

Each elementary conversion decimates the approximating mesh M 1 . The 
30 quality of the approximation therefore deteriorates during the decimation or at best 
remains invariant when the volume remains unchanged, despite the optimization of 
the points. 

For example, if we consider a sphere, it is known that it theoretically requires 
an infinity of triangles to be perfectly modelled. Once the number of triangles is 
35 restricted, the quality of its approximation made to deteriorate. By contrast, in the 
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plane regions of an object, it is possible to decimate until an optimum mesh is 
obtained (thus on a parallelepiped the optimum is achieved when each face is 
described by two triangles). 

It is progressively sought to encode the meshes and therefore iteratively 
5 remove the vertices by a edge merger operation. After the optimization of the 
positions, another process is therefore developed in which a priority is defined on the 
elementary conversions so as to first of all make the conversions that have the least 
possible effect on the model. 

For this purpose, the notion of cost is defined for each elementary operation, 
10 and a priority queue is sorted out according to this cost. 

Thus, initially all the possible conversions are listed and then the cost of each 
of them is computed (without making the corresponding modifications on the mesh). 
The lower cost conversions are made and then the costs of elementary conversions, 
modified by the previous conversions, are recomputed. It may be recalled that the 
15 edge merger eliminates two faces. This means that, as a consequence, the 
conversions pertaining to these faces are withdrawn from the priority queue. This 
creates one or two new conversions for which a computation is also made of the 
energy cost, deduced from the technique described in Appendix 2. This energy cost 
represents the variation in volume between the mesh before conversion and the mesh 
20 after conversion. 

Since the direct computation of volume between two triangulated surfaces is 
complex to compute, the invention makes use of a heuristic method that can be used 
to simplify this problem. 

Let: 

25 - T; be the vertex merging conversion Xj and X,- of M' ; 

X/ the vertex of M' resulting from the merger (for which the position has 
been optimized) f 

F(X i f ) the faces of M' neighboring the vertex X/ after conversion ; 
V M the set of vertices of M belonging to the faces that have been 
30 intersected during the computation of the orientation of the faces during 

the optimization. 

The cost of an elementary conversion T ; is then expressed in the form: 

C(T i (X i ,X j )) = maxd 2 (V M ,F(X i f )) 
This cost therefore corresponds to the maximum distance from the vertices of 
35 the original mesh (beneath the faces neighboring the vertex X { ) to the neighboring 
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faces X, after conversion. Figure 16 gives a 2D view of the expression of the costs in 
the case of a decimation on a discontinuity. As indicated here above, in three 
dimensions, the segments are replaced by triangles. The conversion 181 consists in 
merging the vertices X; and Xj into one vertex X ; '. The vertices 1S2 1 to 182 4 
5 correspond to the set V M . 

It may be noted that, in a meshed plane, the cost of a conversion would be 
zero, since the zones comprising a redundant information would be decimated as a 
priority. 

Exemplary block diagrams are now presented for the implementation of the 
10 different aspects described here above. 
6.3.4.7 Block diagrams 

• Approximation of meshes 

Figure 18 summarizes the approximation algorithm of a mesh M comprising 
n faces by a mesh M' comprising m faces (with of course m<n). 
15 We shall start first of all by copying the source mesh M, or original mesh in 

the variable M', representing the simplified mesh to be determined (191). In this 
mesh M', a search is then made for all the elementary conversions that can be 
performed (192) followed by a computation (193) of the energy costs of these 
conversions after optimization. 
20 The lowest cost conversion among all the possible conversions is then 

selected and carried out (194). This is the decimation step. 

Then, an updating (195) is done of the energy costs of the conversions in the 
neighborhood affected by the conversion made in the step 194. Finally, the new 
conversions induced by the conversion made (194) are added (196) and the 
25 corresponding costs are computed. 

The three steps 194, 195 and 196 are reiterated (197) n - m times until at least 
one of the stopping criteria is reached. " - - 

• Positional optimization 

Figure 20 illustrates the optimizing of the vertex X ; after the merger of the 
30 edge (X i5 X,). 

A preliminary decimation is made of the following initializing parameters, 
partially as a function of the dynamics of the local geometry around the vertex Xj : 

Epsilon is a convergence stopping criterion that is a fraction of the mean 
length of the edges around the vertex X ; (typically 0.001) ; 
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Maxlter is the upper terminal of the number of iterations, arbitrarily fixed 
(typically 200 iterations); 

Steplnit, StepDec and StepMin are respectively the initializing value, the 
rate of decrease and the lower boundary of the step. Steplnit and 
5 StepMin are fractions of the length of the edge to be merged (typically 

0.1 and 0.01), and StepDec is fixed arbitrarily (typically 0.95) ; 
StepSampling is the sampling step on the triangles used for the discrete 
computation of the integral forming the gradient at each iteration. It is 
deduced from the areas of the faces neighboring the vertex X ; 
10 - NbSharpEdges(X) is the number of sharp edges around the vertex X. 

The algorithm starts with a step 201 of positional initialization, which 
consists of the following operations : 

If (NbSharpEdge(X i ) = NbSharpEdge(X j )) 
then 

1 5 initialize in the middle of the edge to be merged 

else 

initialize on the greatest number of sharp edges. 
A computation 202 of the gradient is then done with the discretization 
operations on the triangles and of orientation of the surfaces. Then, the shift 203 is 
20 done. The vertex shifts by the value of the gradient vector multiplied by the current 
step. 

The changes 204 undergone by the step are then checked. If an oscillation is 
detected with the previous gradient vector, the step is multiplied by StepDec. If the 
step is smaller than StepMin, it is reinitialized at StepMin. 
25 Finally, a conversion step 205 is carried out. If the distance between the 

current position and the previous position is smaller than Epsilon, or if the number of 
iterations Maxlter is reached; the "algorithm is stopped (206). Else, the steps 202 to 
205 are repeated (207). 
6.4.3.8 Results 

30 We shall now show, from a simple example (a parallelepiped), the behavior 

of the method of approximation of meshes according to the invention. The usefulness 
of the parallelepiped lies in the fact that it can be used to verify the quality of the 
approximation on the discontinuities (regular sharp edges and corners) and on the 
planes. 
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On an object of this kind, the optimum is reached when each mesh face is 
obtained by two triangles and the overall shape of the object is perfectly preserved. 

Figure 19 illustrates the case of a decimation on a corner 211 of the 
parallelepiped. The initial position was placed on the corner 211 because the number 
5 of sharp edges is the maximum at this vertex (three around the corner 211 and two 
around the regular sharp edge of the object (vertex 212)). The optimization algorithm 
oscillates and then converges around this initial position 213 since the gradient step 
was initialized at a non-zero value and then reduced as soon as the oscillation was 
detected. 

10 Figure 20 shows a decimation on a regular sharp edge. The initial position 

was placed in the middle of the edge 223 since the number of sharp edges is identical 
(two around each vertex 221 and 222). Several positions located on the regular sharp 
edge of the object will comply with volume invariance and the algorithm converges 
after oscillation around the edge on a position that depends on the initial position. 

15 The middle point 224 of the edge therefore constitutes an efficient initialization in 
this case. 

Figure 21 illustrates another situation of decimation on a regular sharp edge 
in which the initial position has been placed on the sharp edge 234 since the number 
of sharp edges are different (two around the vertex 233 located on the edge and none 
20 around the vertex 232 located on the plane face 235 of the object). Several positions 
located on the regular sharp edge of the object will comply with the volume 
invariance and the algorithm converges after oscillation around the edge 234 on a 
position 236 depending on the initial position. 

Figure 22 shows the case of a decimation on a plane 243. The initial position 
25 has been placed in the middle of the edge 245 formed by the vertices 241 and 242 
since the number of sharp edges is identical (none around each vertex 241 and 242). 
Several vertices located in the plane of the object will comply with" the volume 
invariance. The algorithm, after oscillation on either side of the plane, converges on a 
position depending on the initial solution. The middle point 244 of the edge therefore 
30 again constitutes an efficient initialization in this situation. 

The method decimates the initial mesh by using a priority queue. 
Figure 23 shows twelve successive iterations of the algorithm, numbered (b) 
to (m), from the original mesh numbered (a). 
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It is noted that the edges leading to a volume invariance are decimated as a 
priority. The optimum is reached for the vertices requested on this object since each 
face is then meshed by only two triangles. 

It is observed that the edges meeting two corners of the object are not merged 
since this operation would lead to heavy deterioration in the quality of the 
approximation, since the costs and the corresponding elementary conversions are 
greater than the costs of a merger of edges in the plane or along a regular sharp edge. 
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APPENDIX 1 

PROGRESSIVE MESH ENCODING 

Hugues Hoppes has developed a method working by decimation and 
5 optimization of points. The decimation is achieved by a reversible elementary 
conversion that is easily encodable: the edge merger which assembles two adjacent 
vertices into one vertex, eliminates two faces and then optimizes the position of the 
vertex resulting from the merger (Figure 24). 

The author has used a variational approach by which the process is 
10 minimized, the associated energy functional is defined by an error term and a 
regularization term: 

E = E error + E regul 

This error term characterizes the distance between two meshes M and M 
defined locally, and ensures fidelity to the initial data. The regularization term is used 
15 to comply with the topology of a mesh and provides for the unified nature of the 
solution. In order to initialize the problem, once the face to be eliminated has been 
determined, three minimizing operations (starting from three different initial 
positions chosen so as to accurately initialize the optimization process) are achieved 
(Figure 25) in order to choose the solution corresponding to the minimum energy of 
20 the functional. 

Each conversion Tj that eliminates a vertex of the triangulation and modifies 
the mesh M; in M i+1 generates an increase in the total energy that may be computed. 
The energy cost resulting from this decimation is as follows : 

25 / AE Ti =E Mt+1 -E M , (\) 

\ ART; ^ 0 

At the outset, all the possible conversions on the mesh are listed," and then the 
energy cost for each of them is computed. These conversions are stored in a priority 
queue sorted out according to their estimated cost. During the algorithm, the 

30 minimum cost conversion is done and then the cost of the conversions in the vicinity 
affected by this modification are updated (Figure 26). The newly created conversions 
are added and their estimated cost is computed (Figure 27). The algorithm therefore 
acts iteratively and decimates the mesh M of a vertex and of two faces at each 
iteration. This approach of decimation and optimization is well suited to the context 

35 of the progressive encoding since it is enough to record the sequence of the 



31 



decimations in the reverse direction (that of the refinements) to encode a mesh in a 
progressive form. Thus a mesh M may be represented by a mesh with the resolution 
M 0 and by a sequence of refinements {r 0 , r u r 2 , r n }. 

The error used to characterize the distance between the two meshes is the 
Euclidean distance (or quadratic) distance between the point X { to be optimized of 
the mesh M and the mesh M, which has to be defined. The point located at a distance 
D from Xj in R 3 form a sphere with a radius D : 



d 2 (X i; M) = nnn Y) = min \\X % - Y\\ 2 



(2) 



\\Xh = C£ t X?) 1 t 2 ,Xz*! i (3) 



Since the mesh M is a triangulated surface, this generates a measurement of 
the distance between a point X; and a set of triangles. It is observed that according to 
the definition of the L 2 in R 3 that the point of M located at the minimum distance 
20 from Xj is not always unique. The effective choice of this point depends on the order 
of insertion of the points and the nature of the test during the search for the minimum 
(< or <). The metric L 2 applied to the meshes becomes the metric D 2 since the 
reasoning is based on a point-surface distance. 
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APPENDIX 2 

OPTIMIZATION 

The optimization technique chosen to resolve our minimization problem is an 
iterative method known as the adaptive gradient. 

1. Adaptive gradient 

Our goal therefore is to minimize the above-described functional E. Since, in 
every case, the stationarity of E is a necessary condition of optimality, effectively, 
practically all the methods of constraint-free optimization in R N consist in searching 
for a stationary vector X (VE(X) = 0). This problem is equivalent to the resolution of 
the system of non-linear equations: 

~(X) = 0,Vi = l,...,N (1 ) 

It is possible to seek to directly resolve this system. This leads to Newton's 
method. However, this method may be non-convergent if the starting point of the 
iterations is far too distant from X. Furthermore, it assumes that the function is twice 
continually differentiable and requires the computation of the second derivatives at 
each point. 

That is why the most commonly used methods proceed differently : these are 
iterative methods where a sequence of vectors X°, X 1 , ... are generated, converging 
towards a local optimum of E. 

These gradient methods constitute a class of methods that proceed as follows: 
the procedure starts from a point X° and the gradient VE(X°) in X° of the function E 
to be minimized is computed. Since VE(X°) indicates the direction of the greatest 
increase in E, a shift is made by a quantity y 0 in the direction opposite to -the gradient, 
and the following point is defined: 

x - x - 7oX ITv^M! (2) 

The procedure is repeated and generates the points X°, X 1 , X n , 
according to the relationship : 
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(3) 
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In this class, it is necessary to report the determined-step gradient methods in 
which an a priori choice is made of the value of the shift y k . The drawback of this 
procedure is that the convergence may be very slow. We have therefore chosen to 
5 use the adaptive-step gradient method which consists in diminishing (and 
respectively increasing) the step Yk when the error increases (and respectively 
diminishes) while keeping the step control constraint: Yk G [Y m in> Ymaxl- Finally, as a 
stopping criterion, we have taken it as fixed that the difference between two 
successive errors must be below a given value e . 

10 During the minimizing, the positions of the vertices of the mesh will develop 

so as to minimize the volume between the two surfaces. This development requires a 
computation of the gradient at the iteration k (term VE(X k ) of the equation. 
2. Elementary variation of volume 

The process of optimization is therefore seen as a development of the mesh 

15 M so as to minimize v(M,M). The elementary development of the mesh is defined 
by the field of vectors 5M. By carrying out a parametrization of the surface by u and 
v, the vector 5M can be expressed in the form 5M{u, v) . If the surface is locally 
approximated by its tangential plane, the elementary variation of volume is then the 
parallelepiped generated by the development of the surface area element dudv in the 

20 direction .5M(u,v). (Figure 6). 

It is then possible to express the elementary variation of volume induced by 

the area element dudv : \ldM(u,v) dM(u,v)l j 

SV(M, M, 5M{u, v)) = \\ A ^ } ■ 5M(u, v)\ ■ dudv 

The variation is negative when the vector is oriented towards the interior 
25 of the error volume. Figure 7 shows the particular case where the development of the 
mesh is done in the direction of the normal. Let us now consider the fact that : 

fdM(u,v) ~d&(u,v)l 

[~~di A dv- I ' dudv = v ^ da ^ v ) 

30 with n{u,v) being the unit normal at M(u,v) and da(u,v) = II^A^ll- dudv the 
surface area element. 

The variation in distance between M and M is expressed as follows : 
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This expression cannot easily be used because of the absolute value. We shall 
now introduce the function : T](u,v) = . r](M>M,n{u,v)) = l if the normal is oriented 
towards the interior volume and if not -1. We then obtain : 

jn(n, v)8M(u, v) j —r\{u, v)ft(u, v)5M(u, v) 
whence Sd(M,M,5]k) = - J J r](u,v)n{u,v)5til{u,v)da(u,v) 

3. Model of triangular finite elements 

In order to compute the gradient at each vertex of the mesh at the iteration k 
of the optimizing algorithm, the meshing is parametrized by a model of triangular 
finite elements P 1 : 



(4) 



X ; is a vertex of the mesh defined in R 3 and A,;(u,v) is a form function matched with a 
model of triangular finite elements. A form function is defined on the mesh by means 
of barycentric coordinates as described in the document "The method of finite 
elements, basic formulation and linear problems", by O.C. Zienkiewicz and 
R.L. Taylor (AFNOR Technique, 4th edition, 1991). For this purpose, the vertices of 
the triangles (Figure 8) are numbered and a system of coordinates X x , X 2 and X z 
linked by the following relationships is defined : 

x = Ai - x x + A 2 - x 2 + A 3 - 2=3 

y = Ai • yi + A 2 • j/2 + A 3 - yz 
1 = Ai + A 2 + A 3 

We shall now determine the values of the basic functions: - at the vertex 1, Xj 
= 1, = 0 and X 3 = 0. The X l level lines are equidistant lines parallel to the side 2 - 
3 along which X l = 0. It is possible to express the coordinate X x at a point P by a 
ratio of areas: 

_ aire{P2o) 

1 ~ afre(123) ( 5 ) 

, ai + b t ■ x 4- ci ■ y 
We obtain : 1 2A 

X — ° 2 + &2 ' x + ° 2 ' y 
2 ~ 2A 



A 3 



- Q 3 + Q3 - x + c s -y 
2A 



A = aire(123) = 



det 



1 x x 2/1 



1 x 3 

With : ai = x 2 - 2/3 - %3 ■ y2 

h=V2- V3 
Ci = x s - x 2 

The values a 2 , b 2 , Cj, a 3 , b 3 and c 3 are computed by circular permutation of the 
indices 1, 2 and 3. 

4. Surface development 

The development of the mesh can thus be written as follows: 

5d(M,M,SX 1} SX 2t .. .,5X n ) = ~ Xj/ / V^v)n{u,v)8X i X 1 {u,v)da{u : v) 
the partial derivative of the distance for a vertex X ; is deduced therefrom : 

dXi J J K ,„es«p P (A ; ) 
The gradient algorithm for a vertex X ; is written in the form: 



20 where Yi represents an adaptive step with y- t g [y min , y 0 ]- At the outset, the step is 
initialized at y 0 , then multiplied by k (with k < 1) in the presence of an oscillation, 
namely whenfi§|^-, f^-><£KThe behavior of the optimization algorithm may be 
illustrated at 2d (Figures 10 and 1 1), the development of the vertex X; of the model 
M is linked to the sum of the vectors of orientation of the model with respect to the 

25 original curve M, weighted by the form function which is equal to 1 at X £ and 
decreases by 0 along the neighboring segments of X ; . At 3d, the principle is the same 
with a form function defined by the- triangles neighboring the vertex Xj. - 
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AMENDED CLAIMS : 
1. Method for the simplification of a source mesh M formed by a 
plurality of surfaces defined by vertices, faces and orientations of these 
faces, said method implementing a step of decimation by edge merger, 
5 consisting of the association of an edge to be decimated, defined by two 
vertices (11, 12), with a single vertex (13) so as to obtain a simplified mesh 
M\ 

characterized in that the method comprises a pseudo-optimizing step after 
said step of decimation by merger of a edge, positioning the vertex resulting 
10 from said merger as a function of a criterion taking account of a number of 
sharp edges around each of these two vertices forming the edge to be 
merged, so as to reduce the geometrical deviation between said source 
mesh M and said simplified mesh M'. 
^ 2. Method for the simplification of a source mesh according to claim 1 , 

15 characterized in that said step of pseudo-optimization comprises a step of 
enumerating the sharp edges around the two vertices forming the edge to be 
merged and a step of positioning said resulting vertex, in which the following 
two cases are distinguished: 

if the numbers of sharp edges are the same around the two 
20 vertices, the vertex resulting from the merger is placed in the 

middle of the segment linking said vertices (42, 44) ; 
if the numbers of sharp edges are different, the vertex resulting 
from the merger is placed on the vertex with the greatest number 
of sharp edges (43, 45). 
25 3. Method for the simplification of a source mesh according to any of 

the claims 1 and 2, characterized in that it comprises a step for the selection 
of an edge merger to be made among all the edge mergers possible, taking 
account of : 

at least one piece of information representing the curvature 
30 defined locally around the edge considered ; 

at least one piece of information representing the geometrical 
dynamics defined locally. 
4. Method for the simplification of a source mesh according to claim 3, 
characterized in that said step of selection implements a queue of priorities of 
35 edges to be merged as a function of a priority criterion, said information 
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representing the curvature and then a secondary criterion, said information 
representing the geometrical dynamics. 

5. Method for the simplification of a source mesh according to claim 4, 
characterized in that said selection step manages a threshold of curvature, 

5 only the edges with a curvature below said threshold being considered for the 
application of said secondary criterion, 

said threshold being raised when there is no longer any edge having a 
curvature below this threshold. 

6. Method for the simplification of a source mesh according to any of 
10 the claims 1 to 5, characterized in that said information representing the 

geometrical dynamics belongs to the group comprising : 
the length of the edge considered ; 

a mean of the surfaces of the faces neighboring said edge 
considered ; 

15 - a mean of the lengths of the edges adjacent to the vertices 

forming said edge considered ; 

a combination of the lengths of edges and/or surfaces of faces; 

7. Method for the simplification of a source mesh according to any of 
the claims 1 to 6, characterized in that the decimation is interrupted as a 

20 function of one of the criteria belonging to the group comprising: 
a compression rate achieved ; 

a geometrical complexity achieved, expressed by a number of 

vertices or faces ; 

a threshold of curvature achieved. 
25 8. Method for the simplification of a source mesh according to any of 

the claims 1 to 7, characterized in that it constitutes a step of initialization of a 
method of geometrical optimization of a mesh. ~ - - 

9. Method for the geometrical optimization of a source mesh 
comprising a step of initialization implementing the method of simplification 

30 according to any of the claims 1 to 7. 

10. Method for the encoding of a source mesh (M) according to claim 
9, representing a 3D object, delivering a simplified mesh (M') corresponding 
to said source mesh (M), said meshes being defined by a set of vertices, 
edges and/or faces, 



characterized in that it implements a step of simplification according to any of 
the claims 1 to 7, and then a step of minimization of a volume contained 
between said source mesh (M) and said simplified mesh (M'). 

11. Method for the encoding of a source mesh according to claim 10 
5 characterized in that, each of said meshes being defined by the position of 
each of its vertices, said minimizing step provides for the determining of the 
position of the vertices (X-j, X2, X n ) of said simplified mesh (M') 
minimizing the volume V(M, M') between said source mesh and said 
simplified mesh. 

10 12. Method for the encoding of a source mesh according to any of the 

claims 10 and 11, characterized in that said minimizing step implements an 
iterative process progressively optimizing the positions of the vertices of said 
simplified mesh (M'). 

13. Method for the encoding of a source mesh according to claim 12, 
15 characterized in that said iterative process is interrupted when at least one of 
the following stopping criteria is achieved : 
a maximum number of iterations ; 

a difference between two successive shift vectors of the positions 
of the vertices that is below a predetermined threshold (e). 
20 14. Method for the encoding of a source mesh according to any of the 

claims 11 to 13, characterized in that said step of minimization implements 
an adaptive gradient method. 

15. Method for the encoding of a source mesh according to claim 14, 
characterized in that said adaptive gradient method relies on the following 
25 operations: 

- the selection of a vector Xp of R 3n (n > 1 ) of said simplified mesh 
and the computation of the gradient VE (X p ) in X p of the function 
to be minimized E = d(M,M' (X-j,... X n )): 

the determining of the position of X p * of X p of said mesh 
30 according to the relationship defined in the iteration k+1 by : 

k+i k LE(X p k ) 

k varying from 0 to n-1 (with n < N) and yi< being the step of said 
relationship. 
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16. Method for the encoding of a source mesh according to any of the 
claims 10 to 15 characterized in that, at each iteration, an elementary 
variation of said volume corresponding to a vector field SM is determined 
and in that, since the surface is parametrized by u and v so that a vectorSM 

5 is expressed in the form SM (u, v), said elementary variation is likened to the 
parallelepiped generated by the evolution of the surface area element dudv in 
the direction SM (u, v). 

17. Method for the encoding of a source mesh according to any of the 
claims 10 to 16, characterized in that said simplified mesh is parametrized by 

10 means of a model of finite elements. 

18. Method for the encoding of a source mesh according to claim 17, 
characterized in that said finite elements are advantageously obtained by 
means of a refined interpolator. 

1 9. Method for the encoding of a source mesh according to any of the 
15 claims 10 to 18, characterized in that said method implements a progressive 

encoding of said simplified mesh by decimation and local optimization. 

20. Method for the encoding of a source mesh according to any of the 
claims 10 to 19, characterized in that it comprises a step of limitation of the 
deterioration due to an elementary conversion implementing a priority queue 

20 on the elementary conversions. 

21. Method for the encoding of a source mesh according to claim 20, 
characterized in that said step of limitation of the deterioration due to an 
elementary conversion, namely an edge merger, defined by two vertices, 
comprises the steps of : 

25 - computing a cost for each possible elementary conversion ; 

carrying out the lowest cost elementary conversion ; 
recomputing the costs of the elementary "conversions modified by 
the previous elementary conversion ; 

adding the new elementary conversions created and computing 
30 the corresponding costs. 

22. Method for the encoding of a source mesh according to claim 21 , 
wherein the cost of an elementary conversion (Tj) is expressed by: 

C(T i (X i ,X j ))-maxd 2 (V M ,F(X[)) 

with : 
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Tj conversion merging two vertices Xj and Xj of said simplified mesh 
M" ; 

Xj f the vertex of said simplified mesh M' resulting from said 
conversion ; 

5 F(Xj f ) the faces of said simplified mesh M' neighboring the vertex Xj f 

after said conversion ; 

V|\/| set of the vertices of said source mesh M belonging to the faces 
having been intersected during the computation of the orientation of 
the surfaces during said minimization. 
10 23. Application of the method for the encoding of a source mesh 

according to any of the claims 1 to 22 to at least to the following fields : 
virtual reality ; 
scientific simulation ; 
modelling. 
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ABSTRACT OF THE DISCLOSURE 

METHOD FOR THE ENCODING OF A SOURCE MESH WITH OPTIMIZATION OF THE 
POSITION OF A VERTEX RESULTING FROM AN EDGE MERGER, AND 
CORRESPONDING APPLICATIONS 

The invention relates to a method for the simplification of a source 
mesh (M) implementing a step of decimation by edge merger and comprising 
a pseudo-optimizing step after said step of decimation by merger of a edge, 
positioning the vertex resulting from said merger so as to reduce the 
geometrical deviation between said source mesh M and said simplified mesh 
M'. Advantageously, said step of pseudo-optimization consists in 
enumerating the sharp edges around two vertices forming the edge to be 
merged and distinguishing the following two cases : 

- if the numbers of sharp edges are the same around two vertices, 
the vertex resulting from the merger is placed in the middle of the 
segment linking said vertices (42, 44) ; 

- if the numbers of sharp edges are different, the vertex resulting 
from the merger is placed on the vertex with the greatest number 
of sharp edges. 

Figure 4. 
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Legendes des dessins: 

Figure 2 : 

5 curvature (rad) 

Figure 4 : 
Initial position 
Edge to be merged 

10 

Figure 14 : 

Hauteur reportee : transferred height 
Hauteur = height 
Origine = Origin 
15 Cote longueur maxi : Max. length side 
Echantillonage : Sampling 

Figure 15 : 

Relative orientation of the surfaces 
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Figure 16 : 

Vertices to be merged 

Figure 17 : 

25 191 : Copy : M M' 

192 : Search for elementary conversions that can be performed 

193 : Computation of the energy costs of the conversions after optimization. 

194 : Performance of the lowest cost conversion 

195 : Updating of the energy costs of the conversions in the neighborhood affected 
30 by the preceding conversion. 

196 : Addition of new conversions and computation of the corresponding costs. 

Figure 18: 

35 201: Initialization in position 

202: Computation of the gradient 
203: Shift 

204: Changes in the step 
205: Convergence test 
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Figure 23 : 
Iteration 
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Figure 25 : 
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Faces to be eliminated 
Initial positions 

Figure 26: 

Conversion to be updated 
Modified edge 
Merged edge 
Vertex to be optimized 

Figure 27: 
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